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ABSTRACT 



The mass of the Galactic bulge, Mb, is estimated from the tensor virial 
theorem. By including the effects of a barred stellar distribution and figure 
rotation (8f km s~^kpc and by assuming that the bar is oriented at ^ = 20° 
to our line of sight, Mb is found to be as high as 2.8xf0^° M©. This estimate 
is in good agreement with the total mass inferred from the observed optical 
depth to microlensing for stars toward the inner Galaxy. For larger angles, 6 ~ 
40°, or smaller pattern speeds (20 km s~^kpc ~^), Mb is found to be ~ f .7 — 
f.9xf0^° M©, similar to previous estimates for a barred bulge (2xf0^° M©) and 
axisymmetric bulge (f.8xf0^° M©). 



Subject headings: Galaxy: center — Galaxy: kinematics and dynamics — 
Galaxy: fundamental parameters — Galaxy: structure 
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1. INTRODUCTION 

Until recently, the axisymmetric mass model of Kent (1992) was viewed as the standard 
Galactic bulge model; it successfully fits observed kinematics and the near infrared surface 
brightness distribution over the inner few kpc of the bulge. However, analysis of existing and 
new photometric observations show that the bulge may be strongly triaxial. Blitz & Spergel 

(1991) found that the near infrared surface brightness distribution is asymmetric about / = 
0°. Data from the DIRBE experiment on board the Cosmic Background Explorer (COBE) 
confirmed this result (Weiland et al. 1994). The light asymmetry has been interpreted as 
an asymmetry in the stellar distribution. Dwek et al. (1994) fit triaxial models to the light 
distribution and find that the bulge is strongly barred and has its long axis oriented at an 
angle ^ of ~ 10° to 40° to our line of sight. Dynamical modeling of observed gas kinematics 
in the inner Galaxy (Binney et al. 1991) also points to a barred mass distribution with 6 ~ 
16°. 

In addition to its inability to account for the observed light asymmetry, the Kent 

(1992) axisymmetric model appears to significantly underestimate the optical depth to 
microlensing for stars observed toward the inner Galaxy (Paczyhski et al. 1994; Han & 
Gould 1994a; Zhao, Spergel, & Rich 1995). The Kent (1992) model leads to a prediction for 
the optical depth to microlensing which is a factor of 3—4 lower than the observed optical 
depth reported by the Optical Gravitational Lensing Experiment (OGLE) (Udalski et al. 
1994) and the Massive Compact Halo Object (MACHO) experiment (Bennett et al. 1994). 
A barred geometry increases the bulge's optical depth to self-lensing because the optical 
depth scales as the distance from source to lens. A barred geometry naturally places large 
numbers of stars at the near end of the bar (lenses) in front of stars at the far end of the 
bar (lensed sources). Accounting for a barred bulge increases the predicted optical depth 
by a factor of nearly 1.5 compared to the value for the Kent (1992) axisymmetric model 
(Han & Gould 1994a). 

Barred models of the bulge produce predictions for Mb which bracket the estimate for 
an axisymmetric bulge. Zhao (1994) and Han & Gould (1994a) estimate Mb to be 2x10^° 
M© and 1.6x10^° M©, respectively, for barred models of the bulge. This may be compared 
to Mb = 1.8x10^° Mq for the Kent (1992) model (Kent does not state a value for Mb, but 
using his density law and derived mass to light ratio leads directly to this value). In this 
Letter, I demonstrate that a triaxial model of Dwek et al. (1994), but modified to account 
for the excess luminosity in the central 3°, and a predicted value of the bar pattern speed 
(Binney et al. 1991) lead to a bulge mass which may be a factor of 1.4 higher than that of 
Zhao (1994). In what follows, I adopt a value for the distance between the sun and Galactic 
center (Ro) of 8 kpc. 
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2. Mass Estimate 

Given a density model and observational estimates of the kinetic energy and pattern 
speed of a rotating, barred stellar system, an estimate of the system's mass can be made 
from the tensor virial theorem. In the case of the Galactic bulge, sufficient kinematic 
observations exist to estimate the global kinetic energy of the system (see below), recent 
density models have been computed from the observed light distribution (Dwek et al. 
1994), and an estimate of the pattern speed of the bulge potential has been made from 
observations of gas kinematics (Binney et al. I99I). 

For a steady state, rotating, triaxial system, the tensor virial theorem results in the 
following balance between the system's kinetic and potential energy (Binney 1982): 

— Wxx = STj^j; + H^-r^ + iO^{Ix — ly) 

-Wyy = 2Tyy+Tlyy-iO^{h-Iy) (l) 



where Wa^Tn, and Ha [i = x^y^z) are the components of the potential energy tensor, 
ordered motion kinetic energy tensor, and random motion kinetic energy tensor, respectively. 
These equations are for an inertial reference frame with the major axis along the x axis 
and figure rotation with pattern speed io about the z axis. The cu^ term arises because the 
moments of inertia, and ly are not constant in time with respect to the inertial frame. 
This rotation adds an effective kinetic energy along the major axis, and takes an equal 
amount away from the intermediate axis. I have assumed there is no ordered motion in the 
z direction. 

I define the density distribution as p = MBi{x,y , z), where Mb is the system mass 
and is the normalized density distribution. I ignore the disk mass; the effect of this is 
discussed in the next section. Next, I use the following results from Binney & Tremaine 
(1987): 



Ix = i px^dV 
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where V is the volume, r and r' are spatial position vectors, and analogous equations hold 
for y and z. Equation (1) may then be written as, 

-^GM^B^yy = Meiayy^ + ) " Mfi ^ ( z , - ly) (2) 

where wu = Wii/{0.5GMg) and ii = li/Ms- Note that the wu and ii depend 
only on the adopted ^(x^y^z). The kinetic energy terms have been replaced by 
terms involving the system mass and mass weighted mean and random velocities 
(n,, = Jpaf^dV = Meau'; 2Tu = J pvi^'dV = MbIa'). 

To estimate the total mass, I adopt one of the best fitting triaxial distributions to the 
near infrared surface brightness distribution of the bulge from Dwek et al. (1994). This 
is their G2 (Gaussian) model which has principal axis scale lengths of 1487, 584, and 405 
pc (for Ro = 8 kpc). This model is modified inside a major axis radius of 680 pc (415 
pc minor axis radius) to include an axisymmetric density distribution as defined by Kent 
(1992). This hybrid model, similar to one employed by Zhao (1994), accounts for the higher 
observed luminosity in the central bulge which the G2 model does not. The two density 
laws were not smoothly joined, but agree reasonably well at the 680 pc radius. Numerical 
integration yields the values of the wu and ii. 

The pattern speed has been estimated by Binney et al. (1991) to be 81 km s~^ kpc~^ 
from observations of the inner Galaxy gas kinematics. 

The bulge mass is calculated from equation (2) after rotating from the x^y^z frame 
to our line of sight. The angle between the bar major axis and our line of sight is taken 
as ^ = 20°, a value consistent with the G2 model of Dwek et al. (1994). Other models 
which also fit the light distribution well have inclinations up to ^ ~ 40° (Dwek et al. 1994). 
Binney et al. (1991) derive a best fit value of 16°. After rotation, all that remains is to 
adopt a suitable value for the line of sight kinetic energy. Observations at angular distances 
from the Galactic center of 2° — 7° show a large velocity dispersion between about 90 km 
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and 130 km (Rich 1990; Sharpies, Walker, & Cropper 1990; Terndrup, Frogel, & 
Wells 1994; Blum et al. 1994). The mean velocities are generally observed to be small 
in comparison, except at larger longitudes (8° — 10°) where they rise (Minniti et al. 1992; 
Harding & Morrison 1993) to 50—80 km s~^. This increase is compensated for by a decrease 
in dispersion; the total kinetic energy at larger radii appears to remain at ~ (100 km s~^)^. 
I adopt an average of the observed dispersion, (110 km s~^)^, for the line of sight kinetic 
energy. 

The above assumptions lead to a total bulge mass of Mb ~ 2.8x10^° M©. The estimate 
is most sensitive to the adopted angle (6) between the line of sight and bar major axis and 
the pattern speed. The cu^ term is proportional to cos(2^) when rotated from the x,y,z 
frame to the line of sight. If 6 is taken as 40°, the mass estimate falls to 1.9x10^° M©. If 
the pattern speed is 20 km s"^ (and 9 = 20°), Mb = 1.7x10^° M®. 



3. DISCUSSION 



The derived bulge mass is about 1.4 times larger than previous estimates based on a 
variety of techniques. The Kent (1992) dynamical model, which fits nearly all of the existing 
kinematic observations, results in a 1.8x10^° bulge. The non-axisymmetry of the bulge 
coupled with rapid figure rotation allows for a much increased mass. I have ignored the disk 
contribution, which might result in an overestimate of the bulge mass. Applying the tensor 
virial theorem to the Kent (1992) density distribution and using the same observed kinetic 
energy estimate as for the bar model above, I find an axisymmetric bulge mass of 2.3x10^° 
M©. This suggests that ignoring the disk contribution can lead to a somewhat higher mass 
estimate, 20% in the axisymmetric case. The bar model is more compact than the Kent 
model, so the disk contribution may be less in that case. It is also possible that the adopted 
line of sight kinetic energy is too high. This is less likely since the observed kinetic energy 
appears to be ^(100 km s~^)^ out to ~ 1500 pc (major axis radius), as discussed in the 
previous section. 

Dwek et al. (1994) derive a stellar mass to light ratio for their bar model of the bulge 
using assumptions about the initial mass function and evolutionary history for bulge stars. 
Combined with their (luminosity) density distribution this gives a bulge mass of 1.15x10^° 

(Ro = 8 kpc). Zhao (1994) constructed a self consistent dynamical model based on the 
same Dwek et al. (1994) triaxial distribution used in this paper (and similarly modified in 
the central region). Zhao's dynamical model also fits a number of kinematic observations 
for a total bulge mass of 2x10^° M©. Part of the difference between Zhao's model and 
the virial estimate presented in this work is likely due to the fact that Zhao used a lower 
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value of uj (60 km s~^kpc which significantly reduces the cu^ term in equation 2. Han 
& Gould (1994a) also used the virial theorem to estimate the bulge mass (again, using the 
same modified Dwek et al. 1994 density distribution). However, they obtained a lower 
value (1.6x10^° M©) because they did not include the effect of figure rotation (cu = in 
their model). The cu^ term in equation 2 is approximately as large as the kinetic energy 
term if the pattern speed is 81 km s~^kpc Wada et al. (1994) have recently modeled 
the inner Galaxy with a bar potential and run SPH simulations of the gas behavior under 
this potential. They derive a significantly lower pattern speed (~ 20 km s~^) than Binney 
et al. (1991). In this case, the present mass estimate would be reduced to 1.7x10^° M©, a 
value comparable to that found by Han & Gould (1994a) and slightly less than that found 
by Zhao et al. (1995). 

Since the predicted optical depth to microlensing is proportional to the total bulge 
mass (see, for example, Han & Gould 1994a), the mass estimate derived in section 2 would 
result in a predicted optical depth consistent (within the uncertainties) with the observed 
optical depth reported by OGLE (Udalski 1994) and MACHO (Bennett et al. 1994). For 
example, Han & Gould (1994a) report a predicted optical depth of r = 0.95x10"^ for a 
barred bulge. Increasing this value by the ratio of bulge masses (1.75) and adding the 
amount due to disk lensing (~ 0.5x10"^, Paczyhski 1991, Griest et al. 1991) results in r 
= 2.2x10"^. This value may be compared to the observed optical depth, r 3xl0~^, 
reported by Udalski et al. (1994) and Bennett et al. (1994). The observed optical depth 
may be uncertain by as much as 24 % (Han & Gould 1994b), so the prediction from the bar 
model is consistent with the observations. 

Finally, note that the mass estimate made from the observed line of sight kinematics 
may be used to estimate the transverse kinematics. These estimates may be compared to 
existing observations. The longitudinal kinetic energy is predicted to be (153 km s~^)^. 
The proper motion study of Spaenhauer, Jones, & Whitford (1992) results in a longitudinal 
velocity dispersion of 119 km s~^ it 4 km s~^ (Ro = 8 kpc) for stars observed toward Baade's 
window [l,b = 1°, —4°). The predicted value appears rather high even though the observed 
value does not include the effects of mean motion (there was no proper motion standard to 
calibrate the zero point of the measurements). The latitudinal dispersion from Spaenhauer 
et al. (1992), which likely has a small mean motion, is 104 km s~^ib 4 km s~^. This may be 
compared to the global prediction of 101 km s~^ from the tensor virial theorem. 



4. SUMMARY 
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The main objective of this Letter is to demonstrate that the bulge mass may be higher 
than previously thought and that if the bulge is strongly barred, its total mass will depend 
sensitively on the amount of figure rotation (pattern speed) and its orientation to our line 
of sight. 

For the parameters adopted here, the mass of the Galactic bulge is estimated to be 
2.8x10^° M0, about 1.4 times higher than previous bar model estimates. The mass estimate 
is consistent with the mass which would be inferred from the observed optical depth to 
microlensing for stars in the bulge. The mass estimate may be slightly overestimated since 
I have ignored the contribution of the disk in the calculation. A similar calculation by the 
tensor virial theorem for the Kent (1992) model results in a mass estimate which is 20% 
higher than Kent's dynamical estimate. 

The bulge mass derived here depends on factors which are not precisely known, but 
which are observationally constrained. The particular density distribution that I have used 
is only one of a number which fit the observed near infrared surface brightness distribution. 
The mass estimate is most sensitive to the pattern speed and angle between the line of sight 
and bar major axis. The adopted pattern speed is not directly observed through the motion 
of stars, but it results from a dynamical model which fits the observed gas kinematics well. 
If the pattern speed is significantly less than the value adopted here (81 km s~^kpc ~^) or 
the angle between the line of sight and bar major axis is significantly larger than 20°, then 
the mass estimate would be similar to previous estimates for a barred bulge (1.6 — 2.0x10^° 
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